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Abstract 

We consider the singular behaviour of one-loop QCD matrix elements when 
several external partons become simultaneously parallel. We present a new 
factorization formula that describes the singular collinear behaviour directly in 
colour space. The collinear singularities are embodied in process-independent 
splitting matrices that depend on the momenta, flavours, spins and colours 
of the collinear partons. We give the general structure of the infrared and 
ultraviolet divergences of the one-loop splitting matrices. We also present 
explicit one-loop results for the triple collinear splitting, q qQQ, of a quark 
and a quark-antiquark pair of different flavours. The one-loop triple collinear 
splitting is one of the ingredients that can be used to compute the evolution of 
parton distributions at the next-to-next-to-leading order in QCD perturbation 
theory. 
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The high precision of experiments at past (LEP), present (HERA, Tevatron) and future 
(LHC, e^e~ hnear coUiders) particle colhders demands a corresponding precision in theo- 
retical predictions. As for perturbative QCD predictions, this means calculations beyond 
the next-to- leading order (NLO) in the strong coupling as- Recent years have witnessed 
much progress in this field (see Ref. jl] and references therein). In particular, a great deal 
of work has been devoted to study the properties of QCD scattering amplitudes in the 
infrared (soft and collinear) region P] |12j. 

The understanding of the infrared singular behaviour of multiparton QCD amplitudes 
is a prerequisite for the evaluation of infrared-finite cross sections (and, more generally, 
infrared- and coUinear-safe QCD observables) at the next-to-next-to-leading order (NNLO) 
in perturbation theory The information on the infrared properties of the amplitudes 
have also been exploited to compute large (logarithmically enhanced) perturbative terms 
and to resum them to all perturbative orders The investigation of these properties is 
also valuable for improving the physics content of Monte Carlo event generators (see e.g. 
Ref. p^). In addition, the results of these studies prove to be useful beyond the strict QCD 
context, since they can provide hints on the structure of highly symmetric gauge theories 
at infinite orders in the perturbative expansion (see e.g. Ref. [TK]). 

In this paper we consider the collinear limit of multiparton QCD amplitudes at one- 
loop order. We present a general factorization formula, which is valid directly in colour 
space, and discuss some properties of its infrared-divergent contributions. These results 
apply to the multiple collinear limit of an arbitrary number of QCD partons. As an 
example of application beyond the double collinear limit, we present the result of the 
explicit evaluation of a triple collinear configuration of three quarks. Besides its interest 
within the general framework outlined above, our study of the one-loop triple collinear limit 
has specific relevance to the NNLO calculation of the Altarelli-Parisi kernels that control 
the scale evolution of parton densities and fragmentation functions |l6j. This formidable 
NNLO computation is being completed [E| by using traditional methods. Kosower and 
Uwer [18j have proposed to exploit collinear factorization at the amplitude level as an 
alternative method to perform the NNLO calculation of the Altarelli-Parisi kernels. To 
this purpose, the one-loop triple collinear splitting is one of the necessary ingredients. Two 
other ingredients are the tree-level quadruple collinear splitting [TT] and the two-loop double 
collinear splitting. A detailed discussion of the multiple collinear limit and of the results 
presented in this letter will appear in a forthcoming paper |19j . 

We consider a generic scattering process involving final-state QCD partons (massless 
quarks and gluons) with momenta pi,P2, ■ ■ ■ Non-QCD partons (7*, W^, . . .) are always 
understood. The corresponding matrix element is denoted by 

-^a;S:;:^^'^"-(pi,P2,...) , (i) 

where {ci, C2, . . .}, {si, S2, . . .} and {oi, 02, . . .} are respectively colour, spin and flavour 
indices. To take into account the colour and spin structures, we use the notation of Refs. jSl 
1^ . We introduce an orthonormal basis {|ci, C2, . . .) ® S2, • • •)} in colour + spin space, 
in such a way that the matrix element in Eq. (Q) can be written as 

M'^;2;::T'''''-{pi,P2, ■■■) = ((Cl, C2, . . .| ® (Si, S2, ...\) \Ma„a,,...{PuP2, ■■■))■ (2) 

Thus |A^ai,a2,...(Pi?P2, • • •)) a vector in colour + spin space. It is important to specify 
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that the matrix elements A4{pi,p2, ■ ■ ■) that we are considering are physical ones. Their 
external legs are on shell {pf = 0) and have physical spin polarizations. 

The matrix element A4{pi,p2, ■ ■ ■) can be evaluated in QCD perturbation theory as a 
power series expansion (i.e. loop expansion) in the strong coupling as (as = fl's/l^'""))- 
Throughout the paper we are mainly interested in the expansion up to one-loop order. We 
write 



M = {gs 



(3) 



where the overall power q is integer [q = 0,1,2,3,...). In the evaluation of the one- 
loop amplitude Ai^^\ one encounters ultraviolet and infrared singularities that have to 
be properly regularized. We use dimensional regularization in ci = 4 — 2e space-time 
dimensions. The dimensional-regularization scale is denoted by /x. 

The multiple collinear limit of the matrix element in Eq. (Q) is approached when the 
momenta pi, . . . ,Pm of m partons become parallel. This implies that all the particle suben- 
ergies = [Pi+Pj^, with i, j = 1, . . . ,m, are of the same order and vanish simultaneously. 
We thus introduce a pair of back-to-back light-like (P^ = 0,n^ = 0) momenta P'^ and n'^, 
and we write 

iPl + ..-+Pmr = P''+^-^, Sl...m = iPl + ...+PmY , (4) 

2n ■ P 

where Si . m is the total invariant mass of the system of collinear partons. In the collinear 
limit, the vector P" denotes the collinear direction, and we have p^ —>■ ZiP'^, where the 
longitudinal-momentum fractions Zi are 

n ■ Pi 

Zi = ~ (5) 

n ■ P 

and fulfil the constraint Y^'^iZi = 1. To be definite, in the rest of the paper we limit 
ourselves to explicitly considering the collinear limit in the time-like region [sij > 0, 1 > 
Zi > 0). 

In the limit when the m parton momenta pi,...,pm become simultaneously paral- 
lel, the matrix element A4{pi, . . . ,Pm,Pm+i, ■ ■ •) becomes singular. At the tree level, the 
dominant singular behaviour is ^A^^\pl, . . . ,pm,Pm+i, ■ ■ ■) ~ (l/v^)'""^! where s gener- 
ically denotes a two-particle subenergy s^, or a three-particle subenergy Sijk, and so 
forth. At one-loop order, this singular behaviour is simply modified by scaling violation, 
Ai^^\pi, . . . ,Pm-iPm+i, ■ ■ ■) ~ (l/v^)™~^('S/Ai^)~''- The dominant singular behaviour can 
be captured by universal (process-independent) factorization formulae 0-^21 • At the ma- 
trix element level, the factorization formulae are usually presented upon decomposition in 
colour subamplitudes. Nonetheless, collinear factorization is valid directly in colour space 
and we present the factorization formulae in this more general form. 

The colour-space factorization formulae for the multiple collinear limit of the tree-level 
and one-loop amplitudes Ai^^^ and Ai^^^ are: 

l-^Sl...,a„,a„+i,...(Pl, • • . •.•))- 'S'pi°La„ (Pl , • • • , Pm) \M^aL^^,...iP , Pm+1, • • •)) , 

(6) 
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(7) 

These factorization formulae are valid in any number = 4 — 2e of space-time dimensions 
or, equivalently, at any order in the e expansion around d = A. The only approximation 
involved on the right-hand side amounts to neglecting terms that are less singular in the 
multiple collinear limit. Equations (jHI) and ((Tj) relate the original matrix element (on the 
left-hand side) with m + k partons (where k is arbitrary) to a matrix element (on the right- 
hand side) with 1 + k partons. The latter is obtained from the former by replacing the m 
collinear partons with a single parent parton, whose momentum is P and whose flavour a 
is determined by flavour conservation in the splitting process a ^ ai + . . . + am- 

The process dependence of the factorization formulae is entirely embodied in the matrix 
elements. The tree-level and one-loop factors Sp^^^ and Sp^^\ which encode the singular 
behaviour in the multiple collinear limit, are universal (process-independent). They depend 
on the momenta and quantum numbers (flavour, spin, colour) of the m partons that arise 
from the collinear splitting. According to the notation in Eq. Q, Sp is a matrix in 
colour -|-spin space, and we name it the splitting matrix. Making the dependence on the 
colour and spin indices explicit, we have 

gp(c,.^,c™;.,,...,s™)(c.,.„) ^ ((ci,...,C„| ® (Si,...,S^|)5p,^...,^(|c„) ® \Sa)) , (8) 

SO that in Eqs. (0) and ((7j) Sp{pi, . . . ,Pm) acts onto the colour and spin indices ({ci, . . . , c^; 
Si, . . . , Sm}) of the m collinear partons on the left and onto the colour and spin indices 
{{ca, Sa}) of the parent parton on the right. 

The essential difference between Eqs. (jU)) and ((Tj) and the known collinear-factorization 
formulae discussed in the literature jSl El regards the role of the colour quantum number. 
The derivation of the tree-level factorization formula © in colour space is quite straight- 
forward [7]. Its one- loop extension, Eq. ((71), is less straightforward and, in particular, it 
exploits colour coherence of QCD radiation ^^l- The factorization formulae of Refs. [SIE] 
are written in terms of colour subamplitudes and process-independent splitting amplitudes 
denoted by Split(j9i, . . . ,Pm)- Generically speaking, the splitting amplitudes Split can be 
regarded as colour-stripped components of the splitting matrices Sp. For example, in the 
case of the tree-level collinear splitting g ^ Qi + q2, the splitting matrix is 

SpJSf '"^^^^^(pi,P2) = /i^ t^,,, — uiPl) fiP) v{p,) , (9) 

where (c = 1, . . . , — 1) are the SU{Nc) colour matrices in the fundamental represen- 
tation (we use the normalization Tr(t"t^) = Tr5"^, Tr = 1/2), u,v are customary Dirac 
spinors, e is the physical polarization vector of the parent gluon, and the corresponding 
spin indices are understood. The splitting amplitude Splitg^^gj is obtained from Eq. Q by 
simply removing the colour factor t)3^^2. Analogous proportionality relations apply to any 
tree- level and one-loop splitting process a — > Oi -|- 02. When the splitting involves m > 3 
collinear partons, sphtting matrices and splitting amplitudes are not simply proportional. 
Nonetheless, the two formulations of collinear factorization are related by gauge invariance 
and colour algebra, and each formulation has its own practical advantages. 



lMiLa™,a™+i,...(Pl,---,Pm,Pn^+l... 
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Several general properties of the one-loop splitting amplitudes Split*^^''(pi,p2) and, con- 
sequently, of Sp^^\pi,p2) were discussed in Refs. |5ii8jil0j. Here, we present an additional 
general property of Sp^^\pi, . . . ,Pm)- 

When computed in o? = 4 — 2e space-time dimensions, the one- loop amplitude Ai^^^ 
and, hence, Sp^^^ have ultraviolet and infrared divergences that show up as e-poles in the 
expansion around the point e = 0. To make the divergent behaviour explicit, we write the 
one-loop contribution to the splitting matrix as 

Sp^'\p,, ...,Pm)= Sp'^'^^'-ip,, ...,Pm) + ...,pm) , (10) 

where Sp^^'^'^"' behaves as and Sp''^^^^' is finite when e — > 0. We also recall that 
different regularization schemes (RS) can actually be implemented within the dimensional- 
regularization prescription. Two customary RS are the conventional dimensional-regular- 
ization (CDR) scheme (21jj and the dimensional-reduction (DR) scheme (2^. The RS 
dependence of the tree-level and one-loop splitting matrices Sp^^^ and Sp^^^ starts at 0{e^) 
and at 0{l/e), respectively. In Eq. (fTUjl . we define Sp^^^^^' in such a way that it is RS- 
independent at 0{e^), though it is still RS-dependent at higher orders in the e-expansion. 
Therefore, Sp^^^'^"' contains the full RS-dependence of Sp^^^ at ^(l/e^), C(l/e) and 0(e°). 

The divergent part of the one-loop splitting matrix can be evaluated along the lines 
of Refs. |31 EOl 121], and the result can be expressed in terms of a process-independent 
factorization formula. We obtain JTH] 



(4,r)-r(l-2e) 2 \<? -'\ A= J 

1 m 

-E T.-T, (2 - (..r - (.,r) 



~ (t (7. - e^n - (7. - ) - ^ (/?o - epr) 

X , (11) 

where we have used the same notation as in Ref. The colour charge (matrix) of the 
collinear parton with momentum pi is denoted by Tj, and colour conservation implies 
J2i Ti Sp^^^ = Sp^^^ Ta {Ta is the colour charge of the parent parton in the collinear 
splitting). The flavour coefficients 7^ and (3q are 7^ = 75 = "iCpjl and 7^ = /3o/2 = 
(llCyi — 2A^j)/6. The flavour coefficients 7f^ and are RS-dependent [2S1- In particular, 

~CD.R. _ ^GD.R _ ^J^jlg ~D.R _ ~D.R _ ^^^3 aud ^^'^ = fi^'^jl = Ca/Q- 

Note that, so far, we have not specified whether the one-loop amplitudes A^'-^^ and 
splitting matrices Sp^^^ are renormalized or unrenormalized quantities. Since the renormal- 
ization procedure commutes with the collinear limit, the factorization formula (|7j) equally 
applies to both renormalized and unrenormalized quantities. However, the divergent part 
Sp^^^'^"- explicitly given in Eq. refers to the (charge) unrenormalized splitting matrix 
(thus, as is the bare QCD coupling). In the curly bracket of Eq. (fTTj). the contribution 
proportional to Pq — e/^Q^ is of ultraviolet origin, and it would disappear by working at the 
level of renormalized matrix elements and splitting matrices. 
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It is straightforward to check that Eq. (Illj) agrees with the divergent behaviour of the 
known |SJ |HJ ^2] one- loop sphtting amphtudes for the double coUinear limit a ^ ai + a2- 
The triple coUinear splitting process explicitly considered below provides us with a further 
check of the general result in Eq. (|TT|l . 



At the tree level, the splitting amplitudes are known for all possible partonic channels 
in the double |^ ESI ^^"^ triple coUinear limits, and for pure gluonic splitting [g —>■ Ag) 
also in the quadruple coUinear case JT]. In the double coUinear case, the square of the 
splitting amplitudes gives the leading-order Altarelli-Parisi kernels [2Z|. The square of the 
triple coUinear sphtting amplitudes was computed in Refs. [HI EI, and checked in Ref. [TTj. 

At the one-loop level, the splitting amplitudes for all partonic channels in the double 
coUinear limit were computed in Refs. [HI E]] to all orders in e. In the following we consider 
the triple coUinear limit and, more precisely, the coUinear splitting process q — > qi + Q2 + Q3, 
where q and Q denote (massless) quarks of different flavours. 

We first recall the tree-level splitting matrix for the process q q{pi) + Q{P2) + QiPs)'- 

q„(0)(/3l,/32,/33){/3)/ „ „N _ ,,2£V-.C ,C r^N 



X u{p3) Y'v{P2) u{pi) Y u{P) d^y{p2 + P3, n) , 

■Sl23'523 

where j3i{i = 1,2,3) and (3 are respectively the colour indices of the final-state fermions 
and of the parent quark, and 

d^-{k^ n) = -g^'' + 13 

n ■ k 

is the physical polarization tensor of the gluon with momentum k {n^ is the auxiliary 
light-like vector introduced in Eq. (jH)). 

The square of the splitting matrix Sp^^ summed over final-state colours and spins 
and averaged over colours and spins of the parent parton, defines the m-parton split- 
ting function (Pa^...a^), which is a generalization of the customary (i.e. with m = 2) 
Altarelli-Parisi splitting function. Fixing the normalization of the tree-level splitting func- 
tion (Pi°.U„) by 



' \ m— 1 

Si. 



ai ...am ? , (14) 



from Eq. ()12|) we have 

\ giQ2Q3' 2 S23 
where 



■S23'Sl23 Z2 + \ S123 



(15) 



U,,k ^ 2 '''^' '^ ''' + '-^s., , (16) 

Zj^ I Zj Z-f^ I Zj 

which agrees with the result obtained in Refs. [HI Ej . The parameter depends on the 
RS: 5^°-^ = 1 and ^^^^ = 0. 
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To evaluate the one- loop splitting matrix Sp^^^ we use a process- independent method 



[Hinilini- In the case of the collinear splitting process q 
a contribution from two different colour structures: 



Q1Q2Q3 



receives 



4e 



S123 



/33/32^/3i/3'^(Pl'^'2,P3) 



(17) 



The first term in the curly bracket has the same structure as the tree-level contribution in 
Eq. ()12|) . The colour structure of the second term is a new one-loop (quantum) effect. 

The one-loop correction P^^^ to the tree-level splitting function is obtained by simply 
performing the replacement ISp'-'^-'p — > {Sp^^^Y Sp^^^ + {Sp^^^Y Sp^^^ on the right-hand 
side of Eq. (HH). From Eqs. ^ and dHj), we obtain 



qiQ2Q: 



/p(0) 

\^ qiQ2Q3' ^ 2lT 



/ p(l)_(sym.)\ I /p 



(l)(an.)\ 



Note that iP^f}-, n ) is symmetric with respect to the exchange »2 ^ Ps of the Q2 and Q3 
momenta (see Eq. (|15|)). Its one- loop correction has, instead, a symmetric component and 
an antisymmetric component, respectively denoted by (P^^VfJ^ '') and (P^^]^^^^). The anti- 
symmetric component is entirely produced by the one- loop splitting amplitude ^(^1,^2,^3) 
on the right-hand side of Eq. (jl7|) . 

The expression of the splitting amplitude ^(pi,P2,P3) is sufficiently compact to be 
presented explicitly in this letter. We have 



1 , 

J {2'kY 



U{P3 



X df_,u{,q, n) d„p{q + P2+P3,n) 



is2g + iO) (S3g + ZO) 

u{pi)Y{A-4)Yu{P) 



V{P2) 



(g2 + ^0)(tl, + ^0)(s23, + ^0) ' 



where 



tiq = {pi - qf 



{pi + qf 



S23q = {P2+P3 + qf 



(19) 



(20) 



The first and second contributions in the square bracket originate from the one-loop dia- 
grams depicted in Figs. [T^ and[TjD, respectively. 

Note that the expression in Eq. p9|) is valid in any RS. It is also valid in any number 
(i = 4 — 2e of space-time dimensions or, equivalently, to all orders in e. To make the depen- 
dence on e explicit, we have to compute the d-dimensional integral over the loop momentum 
q. This computation requires the evaluation of a set of basic one-loop (scalar and tensor) 
integrals. The corresponding integrands involve, besides the customary Feynman propaga- 
tors 1 /(g^ + iO) , additional propagators of the type l/(n ■ g) , which come from the physical 
polarizations of the virtual gluons (see Eq. (|13|)). Some of these integrals, which resemble 
those encountered in axial-gauge calculations, were evaluated by Kosower and Uwer [TU] in 
the context of their calculation of the one- loop double collinear splitting a — oi + 02. More 
complicated integrals (higher-point functions) of this type are involved in triple collinear 
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Figure 1: One-loop diagrams contributing to the antisymmetric component of the sphtting 
matrix Sp^^Q^Q^. 

sphtting processes. We have computed (to high orders in the e expansion) ah the basic 
one- loop integrals [12] that appear in any triple collinear splitting a — > oi -|- 02 + 03. Using 
these results, we have obtained explicit expressions up to C(e°) of the sphtting amplitude 
A{pi,p2iP3) in Eqs. (fTTj) . (fTI?|) and of the corresponding splitting function {PglQ^Qs'^ 
Eq. (fTH|) . We limit ourselves to presenting the expression of {Pqlo^^Q^}^ since the expression 
of ^(^1,^2,^3) has a very similar structure. 

Including terms up to 0{e^), our result for the antisymmetric component of the one-loop 
splitting function is 



/p(l)_(an.)\ 
\ 91Q2Q3 ' 



r(l + e) T\l 



X 

+ 
+ 
+ 
+ 

X 



S12 
S123 



(47r)"'r(l 
■S123/ 

Z3 

Z2 + Z3 
Sl3 



2e 



4 /-Si23-iO 



4:Nr 



11^ 



Z2 



Z2 + Z^ 
Sl3 



Z3 



Z2 + ^3 



■S13 



6 In 



S123 
■S13 



Z2 



Sl23 " -Sis / V'Sl23- 

S123 S12 ( Z2a- zih 



Z2 + ^3 
S23 
■§123 "~ ■S23 



h\ In 



(0)_ 



S123 

)/iCFTn: 
S23 



S123 



S12 + Si3 — Zi S123 \Si2(l — Zij 

(2gi2 + ^23) a + {si2 - gi3) b 

S12 

ln(^Vnr^)+Li2 (1- 



■S123 / 



S23 



S12 
S123 



Sl3 Z3J \Si23(^2 + ^3), 



-■S123/ \'Si23 

-|- complex conjugate 



Sl3_ 

S123 



+ Li2 1 



■^23 
■5l23 
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(2^3) + 0(e) 



(21) 



where 



S123 / Zi Si3 , ^1(2:1 S23 - Z2 Si3 - Z3 S12) , Zi S23 + Z2 S13 - Z^ Su 



■523 \ -5123 
-5123 / Z2 S23 
S23 V S123 



2(2:2 + 2:3) -512 



2s 



123 



_l_ ^2(^1 -523 — Z2 Si3 + Z3 S12) ^ Zi S23 + Z2 Si3 



Z3 S12 



2(^2 + Z3) S12 

7 



2s 



123 



(22) 
(23) 



It is quite straightforward to check that the divergent part of the one-loop splitting 
function in Eq. ()21|) fully agrees with the result obtained by using the general formula in 
Eq. (fTTjl. Note that the double poles cancel in Eq. (PT|) . so that the most divergent 
terms in {Pq\Q^Q^) are single poles 1/e. These single poles originate from the infrared 
region of the loop integral in Eq. ()19|) . More precisely, they arise when the loop momentum 
q is soft {q 0) but not collinear to any of the external momenta. Note also that, up 
to terms of 0(e°), the full RS dependence of {PqlQ^^q^) is embodied in the corresponding 
dependence (see Eq. (US))) of the tree- level term {P^i^q^qJ on the right-hand side of Eq. (PT|) . 

The absence of additional RS-dependent terms (such as those proportional to 7f^ and 
in Eq. (fTTjl ) at C(e°) is related to the absence of single poles 1/e of ultraviolet and collinear 
origin, as discussed in the second paper of Ref. 



An additional check of our result can be performed by considering the triple collinear 
limit of known one-loop matrix elements, such as the matrix element for the process e+e" 
qqQQ |2H1 121] ■ We have evaluated the one-loop splitting amplitude (fT^ in a spin basis 
of definite ±1 helicities of the quarks and antiquarks. We have then compared the result 
with that obtained by directly performing the corresponding collinear limit of the one-loop 
helicity amplitudes explicitly presented in Ref. |2HI up to 0(e°) (note that the expressions 
in Ref. [2H] refer to the DR scheme). We have found complete agreement. 

The splitting function Pai - -a™, controls the singular behaviour in the multiple collinear 
limit. Moreover, it can have additional singularities (i.e. terms that are not integrable in 
d = 4 dimensions) in some subregions of the collinear phase-space. The tree-level splitting 
function {p'^^q^q^ in Eq. (fT^ is indeed singular when the momenta p2 and pa are parallel 
(i.e. when their relative angle is much smaller than the emission angle of pi) and when 
they are simultaneously soft. The singularity when p2 and are parallel is, instead, 
absent (to any order in the e expansion JH]; and not only at 0{t'^)) in the antisymmetric 
part, {PqlQ^Qs^^ one- loop splitting function. The singularity when p2 and are 

simultaneously soft is still present in {Pqlo'^Q^)^ though its effect vanishes after integration 
over the angles of p2 and p^ because of the antisymmetry with respect to the exchange 
P2 ^ P3- The expression in Eq. (1211) shows that (P^^l,^^n^) has no other singularities in any 
subregions of the phase-space. 

For the sake of brevity, we have limited ourselves, in this letter, to presenting a few 
explicit results for the one-loop triple collinear splitting. These results have mainly an 
illustrative purpose. The method and the tools (in particular, the one-loop integrals) used 
to obtain them are sufficient and can be applied straightforwardly to evaluate the one-loop 
splitting matrix of any splitting process a — ^ ai + 02 + ^s- 
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